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Abstract — Lattice precoding Is an effective strategy for mul- 
tiantenna broadcast. In this paper, we show that approximate 
lattice precoding in multiantenna broadcast is a variant of the 
closest vector problem (CVP) known as ry-CVP. The proximity 
factors of lattice reduction-aided precoding are defined, and 
their bounds are derived, which measure the worst-case loss in 
power efficiency compared to sphere precoding. Unlike decoding 
applications, this analysis does not suffer from the boundary 
effect of a finite constellation, since the underlying lattice in 
multiantenna broadcast is indeed Infinite. 

L Introduction 

Broadcast is referred to as the application where a single 
transmitter sends different messages to many users simul- 
taneously. It may arise, for example, in the downlink of a 
multiuser communication system where the base station wants 
to communicate with the users in the area of coverage. The 
multi-input multi-output (MIMO) technology offers a new 
opportunity for developing efficient broadcast strategies. 

The capacity of a MIMO broadcast channel has been 
determined in yj], where it was shown that the so-called dirty- 
paper coding is instrumental to achieving the capacity. Mul- 
tiple antennas allow to pre-cancel the interuser interference, 
which is known at the transmitter in the broadcast application. 
The lattice method represents a major approach to cancelling 
known interference ||2|l, and specifically, to precoding for 
MIMO broadcast. However, dirty-paper coding suffers from 
high complexity. 

Hochwald et al. yl] formulated precoding as a decoding 

problem at the transmitter Their technique, termed "vector 

perturbation", corresponds to solving the closest lattice vector 

roblem (CVP). It requires the use of the sphere precoder 

, whose average complexity grows quickly with the system 
size. Earher, the idea of precoding via an algorithmic search 
over modulo equivalent points was proposed by Fischer et al. 
for the intersymbyol interference channel |J5|. More recently, 
reference ||6|] considered some practical issues in the imple- 
mentation of vector perturbation. 

To reduce the complexity, lattice reduction (LR) can be 
used, i.e., an approximate solution is found by zero-forcing 
(ZF) or successive interference cancelation (SIC) on a reduced 
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lattice IlZIISI]. Another scheme of approximate lattice precoding 
was proposed in H (we will show that it is actually equivalent 
to LR-aided ZF in (sl]), which was shown to achieve the full 
diversity order 

In contrast to the complexity analysis of sphere decoding 
[lOfl, the complexity of sphere precoding is not available 
in literature. Moreover, the signal-to-noise (SNR) gap be- 
tween sphere precoding and LR-aided precoding has not been 
analyzed, although it has been done for decoding 1 1 III . In 
this paper, we investigate these aspects of lattice precoding 
algorithms. We view the precoding problem as a variant of 
the CVP known as 77-CVP. This view enables us to derive 
the proximity factors for lattice precoding, which measure 
the worst-case loss in power efficiency of LR-aided precoding 
schemes. 

The paper is organized as follows: Section II presents 
the model of MIMO broadcast using lattice precoding, and 
investigates its complexity. In Section III the analysis of the 
proximity factors is given. Section V is a discussion. 

Notation: The transpose, inverse, pseudoinverse of a matrix 
B by B^, B^^, and B^, respectively, and the Euclidean length 
|u|| = -y/(u, u). [.tJ rounds to a closest integer 

II. Lattice Precoding for MIMO Broadcast 

Consider a MIMO broadcast system including one transmit- 
ter, equipped with n antennas, and n receivers, each equipped 
with a single antenna yfl. For convenience, we use the real- 
valued signal model 



y = Hs + n, 



(1) 



where y is the received signal vector at the users, H e M"^" 
is a full-rank channel matrix, s is the transmitted signal, and 
n is the noise vector The entry hij of H indicates the 
channel coefficient between transmit antenna i and user j. s 
is derived from the data vector x = [a;i, ...,x„] . We assume 
that X G [—A/2, A/2]" is taken from the intersection of a 
finite hypercube and an integer lattice. The entries of n are 
i.i.d. Gaussian with variance a'^ each. 



A. Lattice Preliminaries 

An n-dimensional lattice in the m-dimensional Euclidean 
space R™ {n < to) is the set of integer linear combinations 
of n independent vectors bi, . . . , b„ e M™; 



CiB): 




eZ,i = l, 



The matrix B = [bi • • • b„] is a basis of the lattice C(B). In 
matrix form, £(B) = {Bx : x G Z"}. For any point y eM™ 
and any lattice £ (B), the distance of y to the lattice is 
dist(y,B) ~ miiijjgzn Ijy — Bx||. A shortest vector of a 
lattice £ (B) is a non-zero vector in C (B) with the smallest 
Euclidean norm. The length of the shortest vector, often 
referred to as the minimum distance, of C (B) is denoted by 
Ai. 

A lattice has infinitely many bases. In general, every matrix 
B = BU is also a basis, where U is an unimodular matrix, 
i.e., det(U) = ±1 and all elements of U are integers. The aim 
of lattice reduction is to find a good basis for a given lattice. In 
many applications, it is advantageous to have the basis vectors 
as short as possible. The celebrated LLL algorithm is the first 
polynomial (average) time algorithm which finds a vector not 
much longer than the shortest nonzero vector 

Let bi,...,b„ be the Gram-Schmidt vectors corresponding to 
a basis bi,...,b„, where b^ is the projection of b^ orthogonal to 
the vector space generated by bi,...,bi_i. These are the vectors 
found by the Gram-Schmidt algorithm for orthogonalization. 
Gram-Schmidt orthogonalization (GSO) is closely related to 
QR decomposition B = QR. More precisely, one has the 
relations /_* ^ = rij/ri,i and b; = ri^i ■ q^, where q^ is the i*'* 
column of Q. 

A basis B is LLL reduced if 



f^t 



< 1/2 



for 1 < j < i < n, and 



\\Ur>{s-^^l)\\U 



(2) 



(3) 



for 1 < i < n, where 1/4 < (5 < 1 is a factor selected to 
achieve a good quality-complexity tradeoff. 

We now define a variant of the CVP. 

Definition 1 (rj-CVP)): Given a lattice C (B) and a vector 
y e R™, find a vector Bx € £ (B) such that ||y - Bx|| 
<7?dist(y,B). 

B. Sphere Precoding 

In this method, the transmitted signal is given by yD 



s = B(x - Al) = Bx mod C{AB), 



(4) 



where B = H ^, and 1 is an integer vector, chosen to 
minimize the transmission power: 



Note that s e V(£(AB)) (the Voronoi region). The receivers 
apply the modulo function each, obtaining 



y mod A = HB(x - Al) + n mod A 

= (x - Al) + n mod A 
= X + n mod A. 



(6) 



l = argmin|iB(x- Al)p. 

IGZ" 



(5) 



Namely, the data arrive at individual users free of interuser 
interference; the only effect is noise. To solve the CVP (|5]l, the 
sphere precoding algorithm originally proposed for decoding 
purposes was used. 

It is worth pointing out several distinctions between the 
CVP's in decoding and precoding: 

* Decoding gets easier for weaker noise, while noise has 
no impacts on the hardness of lattice precoding. 

• The constellation in decoding is often finite, while the 
lattice in precoding is infinite. Thus, the boundary errors 
in decoding will not an issue in precoding. 

« The received signal in decoding has a Gaussian distri- 
bution centered at a lattice point, while the input to pre- 
coding is roughly uniformly distributed on a fundamental 
parallelepiped. 

For these reasons, sphere precoding incurs more computa- 
tional complexity than sphere decoding at the same dimension 
n. Fincke and Pohst 11211 proposed an algorithm to enumerate 
the lattice points in a sphere, running on an LLL-reduced 
lattice, but their complexity estimate was loose. Kannan's 
algorithm lIlSIl for HKZ reduction can be used to preprocess 
the lattice, giving a CVP algorithm with n""*""'"^ complexity. 
Hanrot and Stehle' improved the CVP complexity analysis to 

7^"/2+o(n) I J4]^ 

On the other hand, Jalden and Ottersten llioll showed that 
the average complexity of sphere decoding is exponential 
with the dimension for any fixed SNR; the constant within 
the exponent, though, does decrease with SNR, meaning 
lower complexity at higher SNR. However, the encouraging 
results for lattice decoding do not extend to precoding. Noise, 
which is crucial to the decreasing complexity of sphere 
decoding, does not even arise in lattice precoding. Since 
the input is largely uniformly distributed in the fundamental 
parallelepiped, the worst-case bound is a sensible measure 
of complexity. Moreover, the paper lllOll assumed a finite 
constellation, rendering the analysis inapplicable to an infinite 
lattice, which is nonetheless the case for precoding problems. 

To conclude, the worst-cast complexity of sphere precoding 
is super-exponential. 

C. Approximate Lattice Precoding 

1) SIC Precoding: To obtain a fast precoder, Windpassinger 
et al. 181] approximated the CVP by using lattice reduction, 
i.e., the closest vector is replaced with Babai's approximations 
||7[]. Let B designate the reduced basis, i.e., B = BU, where 
U is a unimodular matrix. Performing the QR decomposition 
B = QR, where Q has orthogonal columns and R is an upper 
triangular matrix with nonnegative diagonal elements. 



Let u = Q^Bx/j4. An estimate of 1 is then found by the 
SIC procedure: 






,i = n-l,...,l. (7) 



The transmitted signal is given by 

s = Bx - ^Bl. (8) 

At the receivers, the modulo operation is applied, yielding 
y mod A = H(Bx - ABl) + n mod A 



x-^Ul + n mod A 
X + n mod A. 



(9) 



2) ZF P recocting: Let u = Bx/A. An estimate of 1 is found 
by ZF as follows 



1= 



B 



■'H 



A 



The transmitted signal is given by 
s = Bx - ABl 

= BlU-^x-A 



U-^x 



A 



B(U~ix mod A) 



(10) 



The second line of (fTOl i represents the transmission scheme 
in ||8[], while the third line corresponds to the transmission 
scheme in ||9(] . Therefore, the schemes proposed in [jsl |9|] are 
equivalent. To the best of our knowledge, this equivalence is 
not known in Uterature. At the receivers, the modulo operation 
is applied, yielding the same as (|9]l. 

D. Reduction Criteria 

To summarize, the purpose of approximate lattice precoding 
is to find a sub-optimal solution 1 that can reduce the norm 
||s||. Withe lattice reduction, the transmitted vector s falls 
into the fundamental parallelepiped (for ZF) or the rectangle 
spanned by the Gram-Schmidt vectors of the reduced basis 
AB (for SIC). In both cases, the transmission power is 
proportional to the second moment over the respective regions. 
Let V = I dct B| be the fundamental volume of C{B), and V 
be its fundamental parallelepiped. Let ||B|p be the Frobenius 
norm of B, and B be the Gram-Schmidt matrix for B. 
Using a uniform-distribution approximation, the transmission 
powers associated with the approximate lattice precoders are 
respectively given by 



A^ 

PZF = 

Z.t y 



Az A 2 

|x|prfx=^^l|b.|p = ^||Bf 



for ZF, and 



Ps 



IC 



42 " ^ a2 ^ 



4=1 



12 



(11) 



(12) 



for SIC. Therefore, the objective of lattice reduction in this 
application is to minimize the Frobenious norm of B or B. 
However, it is computationally hard to exactly accomplish this 
objective. Thus, the LLL algorithm is often used. 

III. Proximity Factors 

We want to understand the performance of approximate 
lattice precoding. To do this, we compare the transmission 
powers with that of sphere encoding, under the condition that 
they have the same error performance, namely, (0 and (|9) hold 
at the receivers. This is a standard approach to calculating the 
"coding gain'tj. 

The transmission power of sphere precoding is given by 

PsP = ^y'||x||2dx^AV2(V) (13) 

where V denotes the Voronoi region of C{B), and fT^(V) is 
the second moment of V. Then, the SNR gap is asymptotically 
given by 



ll^ll^ for STC 



(14) 



Unfortunately, it is difficult to compute p, and we resort 
to the proximity factors of LLL reduction-aided precoding, 
which measure the worst-case loss in power efficiency relative 
to sphere precoding. More formally, we define the proximity 
factor as 

14^, forZF; 

tW' (15) 

supJ^r, for SIC. 



Fp^ 



sup ■ 



Obviously, p < Fp. This viewpoint implies that the precoding 
problem is jy-CVP: 



|B(x - Al)|l < ri min |iB(x - Al)|l 



(16) 



and consequently, Fp < jf. Babai derived the value of ry ||7|] 
in the case of 6 ~ 3/4. In what follows, we will derive the 
bounds in the general case. Let a = l/((5 — 1/4). 

Lemma 1: If the lattice basis is LLL-reduced, then SIC 
solves 77-CVP for 77 = ?]„ = a'^l'^ l\Ja. — 1. 

Proof: Let B be a LLL reduced basis and B = Bp^ be 
the GSO of the basis B. Given a vector y G M™, we write y 
as a linear combination of the GS vectors y = X]r=i ^i^f Let 
u —'Y^l^xVi^^i be the nearest neighbor of y in £ (B). Let B 
be the integer nearest to /3„ and y'= X]"Ji f^^i + ^b„, and 
V ~0\iri- For n — \, SIC can find the closest vector u. For 
n > 2, we have 



y II = 



/3j||b„||< 



|b„,| 



(17) 



'in practice, power normalization is applied at the transmitter in vector 
perturbation (31l, yet such a scaling factor has no impact on the "coding gain". 



and 



ly-u|i 



\ *=1 



> l^„-Pn|||b„| 

> \Pn-0\\\^n\\ 

= lly-y'll- 



Let w be the lattice point found by ZR Then 

n 

w-y=^/?,b,, (27) 

where |/3j | < 1/2, for 1 <i <n. Let u be the nearest neighbor 
of y in £ (B). We may write 



(18) 



Let w be the estimate of u found by SIC. From ( fTTJ i. we 
obtain 

2 



51 '^i^^ 



(28) 



i=l 



l|y-w|| < 



According to (O, we have 



-Fl|b, 



i=l 



(19) where (/)j e Z. We assume u 7^ w. Let H^^.b/, 



maxi ||(/)jbj||. Then 



|u — wll < n Wdi.h 



fcDfel 



(29) 



1 /q," _ 1 ^ 

1 1 y — w 1 1 < - -1 / 1 1 b„ ! I . 

If pn = 0, then 

||y' - w|| < 7?„_i lly' - ull < 77„_i lly - u|l . 
By (m and (EB, 

||y-w|| = (||y-y'||' + ||y'-wr'''^' 

< (i + r?Li)'/'lly-u|| 

< jynlly-ull. 

If p„ ^ 6*, then 

lly - ull >^^ 



(20) Meanwhile, 



(21) 



where 



u-y = (u-w) + (w-y) 
= (0fe+/3/c) (b/c - m) , 



■"^-^S**'^^'''" 



By (|26l) and |/3J < 1/2, we have 



\^k 



(22) 



Combining this inequality with (|20] i, we obtain 



/a" - 1 

||y-w|| < W -||y-u|| 



l|u-y|l > - 

Combining ( |29] l and ( [30l l, we have 

||u-w|| < n||?ifcbfc|| 



lb, I 



(30) 



vn-l 



< 



a- 1 



< 2n(3V^/2) ||u-y||. (31) 



\/^^T 



|y-u|| 



(23) It is easy to see that 



For a = 2, /y^ = d'^l'^ j \Ja — 1 reduces to Babai's upper 
bound 2"/2. 

Lemma 2: If the lattice basis is LLL-reduced, then ZF 
solves 77-CVP with r/ = ry„ = 1 + 2n (3-ya/2) 

Proof: Let B be a LLL reduced basis. Let Oi be the angle 
between b; and the linear space 5 ([bi,...,bi_i,bi+i,...,b„]) 
spanned by the other n—\ basis vectors. Recall the following 
bound HI 



|y — w|| < ||y — u|| + ||u — w 



< l + 2n(3V^/2)" i|y-u||. (32) 



For a = 2, we have ?;„ = 1 + 2n (9/2)'" ^^^^ 
From the two lemmas, we have the following theorem for 
the proximity factors: 

Theorem 1: According to Lemmas [T] and |2] we have 



sin 9i > 



Since 



we have 



smt/, = mm 



3V"/ 
m bi j 



ie<s lib,; 



(24) 



(25) 



Fp,sic < 



a-l 



and 



Fp.zF < ( 1 + 2n (3^/^/2) 



n-l 



(33) 



(34) 



|m-b,:|| > 



/ 2 



\3y/a 



IbJl , Vm e S. (26) 



These results show that the worst-case loss in power effi- 
ciency of approximate lattice precoders is bounded above by 
a function of the dimension of the lattice alone. 



IV. Discussion 

Our main contribution in this paper was to view the LR- 
aided precoding problem as 77-CVP, compared to the viewpoint 
of bounded distance decoding for LR-aided decoding lll l. 
This viewpoint allowed us to derive the proximity factors, 
which measure the worst-cased bound for approximate lattice 
precoding. Since the underlying lattice is infinite, this analysis 
is rigorous, and it follows that LR-aided precoding also 
achieves full diversity. The derived bounds may not be tight, 
but nonetheless give more insights. Improving the bounds is 
the future work. 
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